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We demonstrate that an arbitrary system of screw dislocations in a smectic- A liquid crystal may be 
consistently treated within harmonic elasticity theory, provided that the angles between dislocations 
are sufficiently small. Using this theory, we calculate the ground state configuration of the TGB a 
phase. We obtain an estimate of the twist-grain-boundary spacing and screw dislocation spacing 
in a boundary in terms of the macroscopic parameters, in reasonable agreement with experimental 
results. 



I. INTRODUCTION 

Condensed matter systems offer a vast stage for the 
intricate interplay between order and disorder. A beauti- 
ful example of such interplay is the twist-grain-boundary 
phase (TGB) of chiral smectics, which is the liquid- 
crystalline analog of the Abrikosov vortex state of type 
II superconductors iQJ^j. Morphologically, the phase 
consists of blocks of pure smectic (which can be either 
smectic- A or smectic-C) separated by parallel, regularly 
spaced twist grain boundaries, where each boundary is 
formed by a periodic array of screw dislocations. The 
direction of dislocation lines rotates by a constant angle 
from one grain boundary to the next. Such a dislocation 
arrangement causes the smectic blocks to rotate about 
the axis perpendicular to the grain boundaries dragging 
the nematic director along. Thus the TGB structure 
combines the properties of smectics and cholesterics: the 
nematic director twists on average as in cholesterics while 
the lamellar structure of a smectic is preserved. In this 
paper only the TGB^ phase will be considered. As sug- 
gested by its name, the smectic blocks in TGB^ are 
smectic- A. 




FIG. 1. Schematic representation of the TGBa phase. 



The analogy between the TGB^ phase and the 
Abrikosov vortex lattice is based on the mathematical 



similarity of the Gibbs free energies for metals in a mag- 
netic field and for chiral smectics jl]-[§| . Their respective 
forms, known as the Ginzburg-Landau free energy and 
the de Gennes free energy, are 
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where ip is a complex order parameter, A the vector po- 
tential, H the magnetic intensity, h a chiral field de- 
termined by molecular structure 0, and n the unit 
director, often decomposed as n = no + 5n. The 



Ginzburg-Landau free energy (1.1) defines two character- 
istic length scales: the order parameter coherence length 
£ = (h 2 /2m*\r\) 1 / 2 and the magnetic field penetration 
depth A = (rn*c 2 g/A-Kne* 2 \r\) 1 / 2 . Their ratio k = A/£, 
called the Ginzburg parameter, controls the phase dia- 
gram as a function of temperature and external magnetic 
field. When the Ginzburg parameter is less than the crit- 
ical value n c = l/\/2, the system is type I, and there is 
a first order transition between a normal metal in a field 
and the Meissner phase with %j) = constant =/= and mag- 
netic field B = 0. When n > n c , the system becomes 
type II, and a new phase intervenes between the normal- 
metal state and the Meissner state. This new phase, the 
Abrikosov flux phase, is characterized by a proliferation 
of linear topological defects of the complex order param- 
eter field ip. The defects are magnetic flux lines, and they 
form a two-dimensional triangular lattice throughout the 
phase. 
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A very similar phenomenon occurs in chiral smectics, 
only the situation becomes somewhat more complicated 
because of anisotropy of the de Gennes free energy (we 
note that the theory of anisotropic superconductors is 
equally complex, see e.g. raj). Instead of a single or- 
der parameter coherence length £, there are two lengths: 
the transverse coherence length £j_ = (C±/\r\y/ 2 and 
the longitudinal coherence length £n = (Cii/jrl) 1 / 2 . But 
since the values of C± and Cii can be made equal by 
rescaling of coordinates, we will not be concerned with 
making distinctions between £j_ and £m and assume £j_ ~ 
£|| ~ £. The chiral smectic analog of the magnetic 
field penetration depth is the twist penetration depth 
A 2 = (K 2 g/2Cq 2 \r\) 1 / 2 , where C » C x « C\\. As in su- 
perconductors, the chiral Ginzburg parameter k 2 = A 2 /£ 
in liquid crystals determines the structure of the phase 
diagram as a function of temperature and chiral coupling 
h. Again, a defect phase, which is now the TGBa phase, 
appears on the phase diagram when n 2 > l/\/2. Linear 
topological defects in the TGB^ phase are screw disloca- 
tions arranged in twist grain boundaries. 

An important problem in the theory of t he d efect 
phases predicted by the Gibbs free energies (1.1) and 
( |1.2[ ) is to relate the lattice parameters of defect lattices 
to the coupling strengths that enter these energies. Thi s 
is rat her complicated because the Gibbs energies (1.1) 
and (1.2) include fouth-order terms that lead to non- 
linear Euler-Lagrange equations. Nevetheless, since the 
publication of Abrikosov's paper over forty years ago 0, 
an extensive body of both experimental and theoretical 
work on the vortex arrangement in the Abrikosov phase 
has been performed, covering all possible ranges of the 
external magnetic field, temperature, and Ginzburg pa- 
rameter ^^j. In contrast with this happy situation in 
superconductors, there is only a small literature that ad- 
dress the liquid-crystal version of the problem, namely 
the determination of the grain-boundary spacing lb and 
the dislocation spacing Id within a grain boundary. In 
part this is due to the fact that rotation of defects in 
TGB^ makes this version of the problem more compli- 
cated. 

In their original paper [jjj Renn and Lubensky consid- 
ered the TGB^ lattice structure near the upper critical 
field h c2 marking the transition from the cholesteric to 
the TGB phase. They employed a model free energy 
with C± = Cii and K\ = K3 = K . Additionally, they 
assumed that the twist pitch is very large in comparison 
to a smectic layer spacing. Under these assumptions, in 
a calculation that parallels Abrikosov's calculation 0, 
they computed the ratio lb/Id for n 2 = 0.80 > 1/V2 and 
various values of K/K 2 . They found that lb /Id increases 
from 0.95 for K/K 2 = to 1.45 for K/K 2 = 10 4 . In other 
words, h /Id is sensitive to the relative values of the Frank 
elastic constants. Our approach, which applies near H c \ 1 
is applicable when the rotation angle between consecu- 
tive smectic blocks is small, predicts a value close to 0.95 
for all values of K/K 2 . 



Experimental determination of the TGB^ lattice pa- 
rameters has proven to be a difficult task. Early exper- 
iments ]To|— |l3|] were crucial in establishing the existence 
of the defect phase of chiral smectics and confirmed that 
the morphology of this phase agrees with the predictions 
of Renn and Lubensky. However, in all these experiments 
lb and Id were estimated rather than measured. Recently 
an extensive structural study of the TGB^ phase occur- 
ring in a series of chiral tolane derivatives ]14| was un- 
dertaken by Navailles and coworkers |ll| . They observed 
an X-ray diffraction pattern in the transverse direction 
to the pitch axis that consisted of discrete Bragg spots, 
which provided a direct way to measure the rotation an- 
gle A0 of smectic blocks. Since the smectic layer spacing 
d is easy to measure in the same experiment, it has now 
become possible to obtain the precise values of Id via 
d/2ld = sin(AO/2). As an additional benefit, the infor- 
mation regarding the number of smectic blocks per pitch 
provided a way to compute the block size from the pitch 
values measured in a different set of experiments based 
on the observation of the Grandjean-Cano steps E3]. The 
ratio lb/Id was found to vary from 0.74 to 1.08. Whereas 
the exact value lb/ Id is up to interpretation, these exper- 
iments are a clear indication that this ratio remains close 
to 1. 

In this paper the equilibrium dislocation arrangement 
in the TGB^ phase is considered in the limit of low-angle 
grain boundaries. Formally, this limit is appropriate close 
to the lower critical field h c i, marking the transition from 
the smectic- A phase to the TGB^ phase. In this case the 
dislocation density is low, and it is possible to neglect 
the interaction between dislocation cores, so the disloca- 
tion core energy associated with destruction of the or- 
der parameter produces an extensive term in the total 
free energy. Thus, in this case, the lattice structure is 
determined completely by the elastic energy cost of the 
distorted smectic structure induced by the presence of 
screw dislocations. In the language of superconductivity, 
this is the London limit of the Ginzburg-Landau equa- 
tions. For our analysis, this means that the results for the 
structure of a single vortex and the systems of parallel 
vortices obtained in the theory of superconductors may 
be directly translated to the langu age of liquid crystals 
using the Gibbs energies ([□]) and Q as a dictionary. 
We will show that there is a simple way to generalize the 
results for parallel dislocations to the case of interact- 
ing twist grain boundaries. In the low-angle limit, the 
rotation angle is a natural small paramater, and to the 
lowest order in this parameter, the interaction energy of 
two twist grain boundaries composed of discrete linear 
defects turns out to be the same as that of grain bound- 
aries with continuous defect distribution. This observa- 
tion allows us to compute the energy cost the dislocation 
arrangement in the TGB^ phase analytically. 

The paper is organized as follows. Section |nj introduces 
the elastic fields suitable for type II smectics-A. We con- 
struct the elastic free energy functional in the harmonic 
approximation. We demonstrate that if we have a solu- 
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tion for the elastic fields uniform in the direction of the 
layer normal at infinity, it is possible to construct more 
general solutions by superimposing shifted and rotated 
copies of the original solution. In Section III we show that 
in the harmonic approximation, systems of pure screw 
dislocations in smectics are distinguished from other dis- 
location systems by a constraint that sets the smectic 
layers to be minimal surfaces. We compute the distortion 
fields and the elastic energy of systems of parallel screw 
dislocations. In Section IV we calculate the elastic energy 
of an individual twist grain boundary and the interaction 
energy of two twist grain boundaries. In Section [y] wc 
obtain the elastic energy density for the TGB^ phase. 
Adding two extensive terms to it, we construct the to- 
tal free energy density. Its minimization yields preferred 
values for the grain boundary spacing and the disloca- 
tion spacing within a grain boundary for the given val- 
ues of material parameters. We find that in the low-angle 
regime the ratio Id /lb is practically constant over a wide 
range of the control parameters, and its value is 0.95. 



II. SUPERPOSITION OF SOLUTIONS FOR 
DISTORTION FIELDS 

The elastic free energy for sme ctics -A directly follows 
from the de Gennes free energy (1.2) by assuming |?/j| 2 
is fixed. Because the smectic density is rapidly varying, 
it is standard to write down the elastic free energy of 
smectics in terms of the layer displacement field u related 



to the phase <!> of the order parameter tp 
decomposition 

</> = $/<7o = n • x - u(x), 
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where no is a fixed unit vector. It is important to note 
that in principle no can be chosen to point in an arbitrary 
direction and the elastic free energy must be invariant 
under rotations of no. In some sense, picking no is like 
fixing a gauge. Dropping the constant terms, choosing 
no to be in the z direction, and denoting n — no as Sn, 
2<7olV'| 2 C|| as B, 2qQ\ip\ 2 C± as D, we have to quadratic 
order in <5n 
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This form of the smectic elastic free energy is known as 
the harmonic approximation. Note that the de Gennes 



free energy and, consequently the elastic free energy (2.2) 
derived from it are invariant with respect to small rota- 
tions but not with respect to arbitrary rotations. There- 
fore, their validity is limited to only those smectic con- 
figurations where the smectic layer normal N and the 
nematic director n do not deviate significantly from n . 



The Euler-Lagrange equations derived from the free 
energy (2J2) are linear, so it is easy to construct new 
solutions by superimposing ones already known. To do 
this, we wil l exploit an underlying symmetry of the full 
theory (1.2). Under an arbitrary rotation R, (f> and n 
transform as 



0'(x)=0(Rx), 
n'(x) = R- 1 n(Rx). 



(2.3) 
(2.4) 



The displacement field u' inherits its transformation 
properties through its definition u = no-x — <fi. In general 
this transformation is nonlinear in the rotation angle 9. 
However, since we are interested in s mall ro tati ons, we 
may expand the transformation laws (2^3), (2^) in the 
rotation angle 9 and keep only the linear terms. 

Consider the distortion fields u and Sn produced by a 
linear source perpendicular to undistorted smectic layers 
at infinity. We will construct rotated solutions by mod- 
ifying simpler solutions: without loss of generality, we 
initially take our source to point along the z-axis. If we 
choose n = z, the distortion fields will be independent 
of z: 



tj){x,y,z) = z- u(x,y), 
Sn(x,y, z) = Sn(x,y). 



(2.5) 
(2.6) 



We may construct a solution which corresponds to a 
superposition of sources that intersect the xy-plane at 
(x a ,y a ) and oriented along = i + 9 a t a , where t Q is 
a unit vector orthogonal to i and 9" < 1. First, con- 
sider what happens to the original solution if the sys- 
tem undergoes a rigid physical rotation about x-axis by 
a small angle 9. The phase function is a scalar, so the 
transformed function is just the original function of the 
transformed coordinates: 

<p'(x, y, z) = <j>{x', y', z') = 0(ar, y - 6z, z + 9y). (2.7) 

In terms of the original elastic field u, the transformed 
phase function is 

4>'{x, y , z) = z + 9y — u(x, y — 9z). (2-8) 

We can read the decomposition ( |2.8| ) in two ways. In the 
first version, it defines the rotated elastic field u'(x,y, z) 
with respect to n' = z + 9y: 

u'(x,y,z)=u(x,y-9z). (2.9) 

The corresponding transformation of Sn is 

Sn'(x,y,z) = Sn(x,y-9z). (2.10) 

If we consider the effect of rotations in Fourier space, we 
note that if 



g(x,y,z) 



d 3 q 
(2tt) 3 



g(q x ,q y ,q z yi* x e^ye^ z , (2.11) 
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then it follows that 



j'(qx,q v , q z ) ~ g(q x , q y ~ Oq z ,q z + 9q v ). 



(2.12) 



The second interpretation of (2.8) is as a definition of 



the rotated field u" with respect to the original no 
z. In this case, the transformed elastic fields acquire an 
extra asymptotic term, linear in y: 



u"(x, y, z) = -6y + u(x, y - 6z), 
Sn"(x,y,z) = 8y + Sn(x,y — Qz). 



(2.13) 
(2.14) 



Even if we have a superposition of sources of different 
orientation, the elastic fields v!' a , Sn'^ produced by indi- 
vidual sources are all defined with respect to the same 
fiducial choice of no = z, so it is consistent to super- 
impose these fields and claim that this superposition of 
fields is the solution for the superposition of sources: 

u sup et(x, y, z) = (2-15) 
J2[ -6 a t a ■ x + u(xjl - xj 1 - 6 a t a z)], 

Oi 

Sn supcl (x 1 y 1 z)= (2-16) 
^[6> Q t a + (5n(x! -x£ - 9 a t a z)]. 

a 

where ii SU per an d <5n super are the total displacement and 
director fields, respectively. The asymptotic parts of the 
distortion fields are absolutely essential in formulating 
global geometric constraints on smectic configurations, 
in particular, for distinguishing different types of disloca- 
tions. When constructing linear superpositions of defects 
we must make sure that the layer normals are unambigu- 
ously defined everywhere. Fortuitously, the terms linear 
in y drop out of the energetics. As a result we are free 
to use either u'(x, y, z) or u"(x, y, z) in our calculations. 
We will exploit this fact in the next section. 



III. INTERACTION ENERGY OF SCREW 
DISLOCATIONS 

A. Screw dislocations vs. edge dislocations 

Dislocations in smectics are linear topological defects 
of the lamellar structure ^ n the presence of dis- 

locations, it becomes impossible to devise a consistent 
global numbering scheme for the smectic layers. In other 
words, the phase field <& cannot be defined as a contin- 
uous field for the whole region where the smectic order 
parameter is defined. If a closed contour surrounds a dis- 
location core and one insists that <J> be continuous, then 
the phase differences along the contour will add up to a 
multiple of 2ir: 



rf$ = 2-7m. 



(3.1) 



However, it is possible to break the space into a number of 
overlapping regions and define $ for each region in such a 
way that local definitions of $ differ in the intersections of 
regions only by a constant. In this case the gradient field 
V$ is globally defined (single- valued and continuous). If 
we choose the same "gauge" no for each region, then the 
field v = — Vj_u defined by V</> = V(&/qo) = no + v 
will be also globally consistent. In terms of v, the above 
integral can be rewritten as 



v • dl = nd. 



(3.2) 



It is convenient to represent a dislocation line of strength 
n and position x = x(^) by a singular dislocation density 



b(x) 



dl 



(3.3) 



With its help, the integral relation (3.2) can be turned 
into a differential form: 



Vx v = b(x). 



(3.4) 



For a system of dislocations, the density is constructed 
by s uperimposing densities of individual lines in the form 
( |3.3|) . Clearly, this procedure guarantees that all integrals 
in the form (3.1) or (3.2) will have the right value. In the 
most general case, distributions of dislocations can be 
continuous. In principle, any function b(x) can be con- 
sidered to be a dislocation density as long as it satisfies 
a conservation law V- b = 0, which says that dislocation 
lines cannot end inside the system. 



The integral co nstr aints (3.1) and (3.2) or their differ- 
ential equivalent (3.4), define dislocations topologically, 
but by their nature they are incapable of determining 
how dislocations are actually arranged in physical sys- 
tems. To understand that, we need to look at disloca- 
tions as physical objects that have energy. There are two 
contributions to the dislocation energy. The core energy 
F'corc arises from the destruction of the order parameter 
in the core region. It is proportional to the total length 
of dislocation lines in the system, but does not depend 
on the details of the smectic configuration outside the 
core. Still, for topological reasons, the smectic structure 
outside the core is necessarily distorted, and these distor- 
tions give rise to the elastic energy F e \ of the dislocation. 

The structure of the distortion fields induced by a dis- 
location depends significantly on the orientation of the 
dislocation core with respect to the smectic layers. This 
orientation reflects the physical nature of the dislocation. 
We may consider the Volterra construction for two classes 
of defects ||. One way to produce a dislocation in a 
smectic is to force the layers to wind in a helical fashion, 
which causes the core to be perpendicular to the smectic 
layers at infinity. A physically different procedure is to 
remove half a layer, which makes the core parallel to the 
layers. Dislocations of the first type are known as screw 
dislocations and of the second type as edge dislocations. 
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Now let us turn to calculation of the elastic energy of 
a system of screw dislocations. The program is to con- 
struct an appropriate dislocation density b(x) and then 
minimize the sm ectic elastic free energy while requiring 
compliance with (3.4). If we restrict our consideration to 
systems of nearly parallel screw dislocations with sepa- 
rations large compared to the preferred layer spacing d, 
then the criteria for the validity of the harmonic approx- 
imation are met and we can use the elastic free energy 

(U). 

First, consider a single screw dislocation at the origin. 
The corresponding dislocation density is 



bi(x) = z dS(x)S(y). 



(3.5) 



If we make the natural choice no = z, the distortion fields 
v and (5n are cylindrically symmetric. It turns out that 
t he cylindrical symmetry and the topological condition 
( |3.4|) completely determine the v-field. 

An immediate consequence of the cylindrical symmetry 
of v and Sn is that their derivatives in the z-direction 
vanish: d z ~v = 0, d z 5n = . T he x- and y-components of 
the topological condition (3.4), imply that d v v z = and 
d x v z — 0. It then follows that v z is constant, set to zero 
by the boundary conditions at infinity: 



= 0. 



(3.6) 



Therefore, the layer displacement field u is independent 
of z, and the distortion fields induced by a single screw 
dislocation are of the type considered in the previous sec- 
tion. Thus, we may construct the distortion fields for the 
entire system of nearly para llel sc rew d isloca tions by the 
superposition procedure of ( 2.15 ) and ( 2.16| ) . 

There is another constraint on the distortion fields im- 
posed by the cylindrical symmetry. Since d z ~v and d z Sn 
vanish, the compression term and the bend term drop out 
from the elastic free energy (2.2). The remaining terms 
involve only the radial and azimuthal components of the 
distortion fields. Moreover, there is no term that mixes 
different components. Observe that the topological con- 
dition ( ft.4| ) fixes the azimuthal component of v. Among 
all configurations of the distortion fields with fixed az- 



imuthal components, the configuration with v p 



0, 



5n p = has the lowest energy. This means that the 
distortion fields induced by a single screw dislocation are 
not only cylindrically symmetric, but also have a single 
azimuthal component. A vector field of such structure is 
solenoidal, and we can conclude that 



V- v = 
V- <5n 



0, 

= 0. 



(3.7) 
(3.8) 



The last relation eliminates the splay contribution to the 
elastic free energy. Although (|3.6|), (|3.7[), and ( |3.8| ) are 
obtained for a single screw dislocation, the superposition 
procedure guarantees that they remain valid even for the 
distortion fields produced by arbitrary systems of screw 



dislocations as long as the harmonic approximation is 
applicable. 



A simple argument shows that (3.7) cannot hold in 
the presence of edge dislocations. Consider a single edge 
dislocation oriented in the y-direction. As in the case 
of a screw dislocation, the translational symmetry of the 
distortion fields in th e co re direction together with the 
topological condition (3.4) makes the v-field orthogonal 
to the core direction, so v y = 0. Then V ■ v reduces to 
d x v x + d z v z . It is easy to check that d x v x and d z v z have 
the same sign everywhere, so V • v never vanishes. The 
presence of additional edge dislocations does not modify 
this conclusion. 

Thus we see that in the harmonic approximation, the 
constraint (3/7) provides a very clear distinction between 
systems of pure scre w dis locatio ns a nd other dislocation 
systems. Note that (3J3) and (3/?) imply that smectic 
layers in the presence of screw dislocations are minimal 
surfaces, which are defined as surfaces of zero mean cur- 
vature. If N is a field of unit normals to a family of sur- 
faces 4>(x,y,z) = const filling the space, then the mean 
curvature of the surfaces in the family is proportional 
to V • N |l8) , the multiplicative constant chosen by con- 
vention. The field of unit normals to the smectic layers 
is 



N = 



V(j> z + v 

= VI + 2v z + 



(3.9) 



Since v z vanishes, |V(/>| — 1 + 0(v 2 ). So in the harmonic 
approximation the denominator is unity, and 



V- N = V- 



0. 



(3.10) 



B. Elastic energy of screw dislocations 

Following Q , we can obtain an explicit solution for the 
distortion fields induced a single scr ew d isloca tion . The 
v-field is deter mine d by equations (3.4) and (3.7) with 
b(x) given by (3J)). At infinity the smectic layers are 
undistorted, so v — > as p — ► oo . To emphasize tha t the 
problem is in fact two-dimensional, we rewrite (3.4) and 
© as 



Vj_ x v = b(x), 
Vl-v =0. 



(3.11) 
(3.12) 



For our purposes, it is convenient to solve these equa- 
tions in Fourier space. A general solution to (3.11) can 
be decomposed into longitudinal and transverse compo- 
nents with respect to the direction of cjx- Equation (3.12) 
eliminates the longitudinal component, so we have 



v(qj_) = i 



.qix b(qj_) 



qi 



(3.13) 
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where in the case of a single screw dislocation b(q) = 
zdS(q z ). The nematic director tilt field Sn can be ob- 
tained from the Euler-Lagrange equations derived from 
the elastic free energy fl2.2j) . Earlier we discovered that 
the compression, bend and splay terms drop out. Then 
the elastic free energy reduces to 

F = \ J d 3 x [D{X7 ±U + Sn) 2 + K 2 {X7 ± x Sn) 2 ] 
d 3 q 



(2tt) 



D|v(qi_) - *n(<U-)l' 



Xalqxx (5n(qL)| 2 



(3.14) 



Variation of ( 3.14 ) with respect to u and Sn leads to the 
following Euler-Lagrange equations: 



Vl- (v-<5n) = 0, 
V|5n= ^(<5n-v), 



(3.15) 
(3.16) 



where A = A2 = (K2/D) 1 ' 2 is the twist penetration 
depth. The first equation is automatically satisfied. The 
solution of the second equation in Fourier space is 



fa(qi) = iTW v(qi) ' 



(3.17) 



In real space, the solutions for the distortion fields in 
polar coordinates p and <fi are 



V= 2^ 



v-<5n= — JC 1 {p/X)e 4 
In \ 



(3.18) 
(3.19) 



where K,\ is the modified Besse l func tion o f orde r 1. Sub- 
stituting the distor tion fields ( 3.13| ) and ( 3.17 ) into the 
elastic free energy (3.14), we find that the elastic energy 
cost per unit length of a single screw dislocation is 



R 



(i) 



1 



L 



d 2 q 



1 



Dd 2 

2 Ua J (2ir) 2 q2+l/A2- 



(3.20) 



Since the smectic order exists only outside the dislocation 
core, the integration region has to be restricted to the 
disk |qj_| < 1/a, where a is the core radius. Then the 
integral gives 



F 



(i) 



Dd 2 , f\ 2 . 



(3.21) 



Note that in the extreme type I systems A — > 0, so there 
is no elastic contribution to the screw dislocation energy 
in the harmonic approximation p| fi"7||. 



It is straightforward to generalize (3.20) to an arbitrary 
density of screw dislocations parallel to the z-axis: 



Fel l n f d 2 q 1 . |S 

- = 2 D J (2^qITW |5(qi)l 



(3.22) 



where the scalar areal dislocation density &(qj_) is de- 
fined via b(q) = zb(q±)S(q z ) = 
zd8(q z ) J2 a exp{iq x x a } ev^>{iq y y a }. 



IV. INTERACTION OF TWIST GRAIN 
BOUNDARIES 

A twist grain boundary separates two smectic domains 
with layer normals that, while pointing in different direc- 
tions, remain perpendicular to some axis. Physically, it 
can be implemented as an array of equidistant parallel 
screw dislocations. In this case topology imposes a con- 
straint on the layer rotation angle AO. Rewriting the 
topological condition ( |3.1| ) in terms of the layer normal 
and the local layer spacing c?(x) Mf, 



N 



(4.1) 



If we consider a rectangular integration path in the xy- 
plane which surrounds one n = 1 defect then, with d 
fixed, we find that 



l j [SN + - JN_] = 1, 



(4.2) 



where 5N is the projection of N onto the xy-plane, d 
is the equilibrium layer spacing and l c i is the dislocation 
spacing. This change corresponds to the rotation of smec- 
tic layers by AO = 2sin _1 (d/2^), which becomes d/ld 
in the low- angle limit. It should be emphasized that al- 
though the rotation angle of the smectic layers is dictated 
by topology, topology in no way requires a specific ori- 
entation of defects in the grain boundary with respect 
to the smectic layers at infinity. Rather, this orientation 
is determined energetically. Leaving a detailed compar- 
ison of energetics of various dislocation systems outside 
the scope of this paper, here we will assume that pure 
screw dislocation systems are energetically preferrable to 
systems of screw-edge dislocations with similar geome- 
try and consider only those configurations of defects and 
smectic layers that correspond to pure screw dislocation 
systems. The formalism developed in the preceding sec- 
tion applies specifically to this kind of system. 

The above result regarding the rotation angle of the 
smectic layers can be also obtained by considering the 
distortion fields. If the plane of the boundary is the yz- 
plane and the dislocations are parallel to the z-axis, then 
the dislocation source has the following form: 

00 

bt g b(a:,2/, z) = zd6(x) 2J S(y - nld), (4.3) 

n— — 00 

where Id is the spacing between defects along the y-axis. 
The layer tilt v induced by the source (4.3) can obtained 



G 



by su perim posing the contributions of individual disloca- 
tions (|3.18|) : 



d ^ y - nl d 



2tt ^ + (y - nl d V 

n= — oo va °' 



2-7T ^-^ + (y — nld) 2 

n— — oo w ' 



(4.4) 



(4.5) 



The sums can be computed explicitly with the help of 
the Poisson summation formula |19fl: 



sin 2wy/ld 



2ld cosh2-Kx /Id — cos27rt//Zd 
d sinh27rx/Zd 
2ld cosh27rx/ U — cos27ry// fl 



The limiting form of v x and v y for large \x\ is 



v(x -> ±00) = ±^-y, 



(4.6) 
(4.7) 

(4.8) 



which shows that the smectic layers undergo a rotation by 
d/ld about the cc-axis as they cross the dislocation array 
at x = 0. While the director relaxes to the layer normal 
in a distance of the order of the twist penetration depth 
A, we see that the smectic layers relax to the undistorted 
asymptotic configuration within a distance Id of the grain 
boundary. 

The intensive energetic characteristic of a twist grain 
boundary is the elastic energy per unit area. It can be 
computed by the following limit ing procedure. Instead 
of an infinite dislocation array (4.3), consider an array 
of N screw dislocations. The array extension in the y- 
direction is Nld- Its elastic energy F e \(N) is given by 
(3.22). Then the elastic energy per unit interval of the 



y-axis of a complete twist grain boundary can be taken 
as a limit of the ratio of the energy of a finite array to its 
extension in the y-direction: 



F 



(i) 



t g b 



Id 



Urn 

N-*a 



F cl (N) 
Nld ■ 



(4.9) 



To implement this limit, we should consider the square 
of the amplitude of the areal dislocation density 



N N 

\b\ 2 = d 2 J2Yl e lq y ld( - n - n ' 



(4.10) 



As N — * 00, note that \b\ 2 /N — > J2™=-oo exp{iq v ld,n} 
and thus the energy per unit area of a twist grain bound- 
ary is 



Dd 2 
~2ld~ 



E 

n— — 00 



d 2 q 



^iq y l d n 



(2n) 2 q 2 



1/A 2 



(4.11) 



Note that F tg b/A may be broken into an extensive part 
and the interaction part. This corresponds to a separa- 
tion of the 11 — term from the rest of the sum: 



- tgb 



1 F, 



(i) 



In 



A ' 



(4.12) 



The interaction part F^ can be written as a sum over 
contributions of dislocation pairs at distances Id, 21 d, ■ ■ ■ : 



Dd 



d 2 q 
(2^ql 



iq y l d r 



1/A 2 



Dd f 00 dq y cos q y ldn 



2l d ^ 



— 2lT 4^ + 1/A 2 
D<P^ r (l d n\ 

n—l v y 



(4.13) 



where /Co is the modified Bessel function of order 0. As 
expected, this result closely resembles the interaction en- 
ergy of parallel vortices in the London limit [^o|,|l| : 



^Abr 



1 $ 



71 £ 8tt 



n?£I>(f) 



(4.14) 



% J>1 



where $0 = 2nhc/e* is a quantum of the magnetic flux, 
til is the vortex areal density, A is the magnetic field 
penetration depth. 

Our next goal is to compute the interaction energy 
of twist grain boundaries. Before we consider the grain 
boundary system in the TGB^ phase, let us limit our 
consideration to systems of finite number of parallel low- 
angle grain boundaries that are suf ficie ntly separated 
so that the harmonic approximation (2.2) is applicable. 
Again, topology completely determines the relative ori- 
entations of the smectic layers in different smectic blocks. 
To ensure that the defects are pure screw dislocations, we 
require that the defects in adjacent grain boundaries be 
rotated by AG = d/ld and, in addition, that the outer- 
most smectic blocks be rotated by the same angle, but 
in the opposite directions with respect to the defects in 
the middle of the system. In comparison with the cal- 
culation for a single grain boundary, we are dealing now 
with more general defect systems where not all defects 
are parallel to each other. Our formalism easily handles 
this situation. 

Consider a system that contains only two twist grain 
boundaries separated by a distance I from each other. 
The angle between directions of the dislocation lines in 
these two boundaries is 9 — AO = d/ld- We may im- 
plement this by taking one grain boundary at x = —1/2, 
rotated by —9/2 and the other boundary at x = 1/2, ro- 
tated by 9/2 where 9 = d/ld- The dislocation density for 
this complexion is the sum of the contributions from the 



two grain boundaries, b = b 



(i) 
t R b 



b[ g b> where 



h^ h {q x ,q y ,q z ) = [z-(8/2)y] 2nde-^ 1 / 2 x 



S[q z - (9/2)q y ] 



(4.15) 



71— — OO 
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b[l(q x ,q y ,q z ) = [z + {fi/2)y] 2vde^ 1 ' 2 x 



6[q z + (0/2)q y ] c 



i(q y —9q z /2)nl d 



(4.16) 



We may transform the sums in the above expression into 
a sum of delta functions via the Poisson summation for- 
mula: 



b&fa*, 9„ 9.) = [5 - (0/2)$] <5[g 2 - (9/2)q y \ x 



2?rd 



-iqxl/2 



m— — oo 



,/2- 



2irm 



bS(&. 9». &) = ^ + (*/ 2 )fl 5 fe + (*/%»] x 

27rm 

— < - > ,^ a - tin. r> - 



J2- 



(4.17) 



(4.18) 



By virtue of linear superposition, the energy of this dis- 
location density will be the sum of three terms. Two of 
these terms are simply the self-energies of the two indi- 
vidual grain boundaries which we have calculated above. 
The interaction energy comes from the cross term which 
is of the form: 



d 3 q h (i),j 
(2tt) 3 *s b 



(q)My(q)&£ J '(-q), (4.19) 



where M y (q) is the general interaction kernel that ac- 
counts for bot h scr ew an d edg e dislocations ||. Upon 
inspection of (4.17) and (4. IS), we note that: 1) Since 
(d/ld) = 9 the y-components of the dislocation densities 
may be ignored since we are only working to quadratic 
order in 9 = A9, and 2) in the product ^gb^tgl on ^ tne 
rn = terms in (4.17) contribute because of the delta- 



function constraints. As a result, only the q z — and 



q y = modes of and contribute. Moreover, 
since the a, 



(2, 



= modes of the layer normals in fl4.6| ) 
and (4.7) relax to their asymptotic values immediately, 
we see that we have constructed a globally consistent 
laye r stru cture between the boundaries. Thus we may 
use (3.22) to evaluate the interaction energy: 



teb 



(2, Dd 2 1 



2ld Id 

Dd 2 X 

2ld Id 



dq x e 



iq m l 



-iq x l 



2tt 



1/A 2 



(4.20) 



Note that since only the q y = modes contribute, the in- 
teraction energy is independent of arbitrary phason shifts 
of the grain boundaries along y. In the harmonic approx- 
imation, the interaction energy of grain boundaries in a 
system of several twist grain boundaries breaks down into 
the sum of contributions of individual grain boundary 
pairs. Since the h armon ic theory admits linear superpo- 
sition, the result (4.20) can be applied to any angle of 
rotation at any separation I as long as the grain bound- 
aries are composed purely of screw defects. 



Another way to look at this result is to consider the 
behavior of the director between the grain boundaries. 
We may expand the director and displacement fields in 
Fourier modes at all the reciprocal lattice vectors G. In 
the harmonic approximation, director and displacement 
fields from different sources add linearly. Thus if no pairs 
of grain boundaries have dislocation axes parallel or an- 
tiparallel and dislocations in each boundary are straight, 
each distortion (displacement u G or £n G ) for a given re- 
ciprocal lattice vector G / arises from a unique grain 
boundary. Thus finite reciprocal lattice vector distortions 
from different grain boundaries do not interact. Each 
grain boundary, however, produces a G = director dis- 
tortion Suq, which is sensitive to the presence of other 
grain boundaries. The origin of interactions between 
grain boundaries is thus 5no, and we can calculate these 
interactions by applying appropriate boundary condition 
to 5no. For an isolated grain boundary, <5no = (0, 6n y , 0) 
reaches constant asymptotic values of (O,±0). If there is 
more than one grain boundary, Stiq has to rotate through 
the angles determined by the dislocation complexion in a 
shorter distance and at greater energy cost. Consider two 
walls with dislocation separation Id located at x = ±1/2. 
If we consider (3.15) and ( p. 16 ) in this situation we see 
that • v = • <5n = since both v and Sri only have 
components along y but only depend on x. Thus ( 3. 16 ) 
becomes: 



(4.21) 



We have three regions to consider. For x < —1/2 we have 



tin,. 



V_i_u = 0y, 

-9 + Ae(*+*/ 2 )/\ 



while between x = —1/2 and x = 1/2 

Viu = 0, 
sinh(a;/A) 



5Uy ~ B "sinh(Z/2A) : 



and for x > 1/2 



Vl U : 



-0y, 



5n y = 9-Ce- ( ~ x - l l 2 ^ x . 



(4.22) 
(4.23) 



(4.24) 
(4.25) 



(4.26) 
(4.27) 



Continuity of the director forces C = A = 9 — B. Insert- 
ing these solutions into into the energy (3.14) we have 



FQ) 
A 



DX [(B - Of + B 2 coth(Z/2A)] 



(4.28) 



Minimizing over the free parameter B and using 9 = d/ld, 
we find 

F(l) _ DXd 2 coth(Z/2A) 



l\ l + coth(Z/2A) 
DXd 2 r 



i 



-l/X 



(4.29) 
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The energy of interaction is simply 
AF F(l) - F(oo) Dd 2 



A 



HI 



(4.30) 



which agrees, with our previous result ( 4.20 ). This shows 
that the energy of interaction comes from the "confine- 
ment" of the director - it is forced to twist from —9/2 to 
9/2 in a length on the order of a few I. 



V. TGB A phase 

In the preceding section we showed that in a sys- 
tem composed of finite number of low-angle twist grain 
boundaries the energy of the dislocation interaction 
within a grain boundary and the energy of the inter- 
boundary interaction decay exponentially with distance. 
It suggests that the dislocation arrangement in the TGB^ 
phase could be treated as well within the same computa- 
tional framework, even though the angles between the 
directions of dislocations in the entire system are not 
restricted to only small ones. Indeed, when the grain 
boundaries are low angle, the dislocations that cannot be 
described as nearly parallel are separated by the distance 
of many grain sizes If,. It is reasonable to hope that we 
can find a regime where lb is sufficiently large so that the 
interaction part of the TGB^ elastic free energy density 
is dominated by the interaction of the dislocations in a 
few nearby grain boundaries. In this case our formalizm 
would be reliable. 

Provided that that we are in the correct regime, the 
preceding section supplies us with all the essential ingre- 
dients to compute the interaction energy of dislocations 
arranged in the TGB^ structure. To construct an appro- 
priate dislocation source, we need to combine the sources 
for individual twist grain boundaries at positions 0, ±Zf,, 
±2lb, . . . along the pitch axis. The grain boundaries at 
different positions are distinguished by the direction of 
defects and, in addition, might be arbitrarily shifted in 
the direction perpendicular to the pitch axis. We may, 
through superposition, use the results from the last sec- 
tion to find: 



fint 
/TGB 



Dd 2 

2lbld 
A 



OO 



71=1 



A 



(5.1) 



We can understand the interaction term by again con- 
sidering the confinement energy of the director. In the 
case of the full TGB^ structure, each "cell" between the 
grain boundaries must be identical. Thus we may use 
(4.24) with B = 9/2 to calculate the energy for each cell. 
We find complete agreement with the interaction term 
above. 




0.05 0.1 0.15 0.2 

smectic | TGB-A a = 2(hd - E)/Dd 2 -5- cholesteric 

FIG. 2. Dependence of the dislocation spacing within a 
grain boundary Id and the grain size h on the control pa- 
rameter a = 2(hd — E)/Dd 2 . At some large value of a there 
is a transition to the cholesteric phase (at hcz) while a — 
corresponds to the smectic-TGBA transition. 



In addition to the elastic energy of interacting screw 
dislocations, the total free energy of the TGB^ structure 
includes two extensive terms that depend only on the 
dislocation density 1 /{Ibid) and not on the details of a 
particular dislocation arrangement. These terms are the 
extensive part of the screw dislocation energy density and 
the chiral energy term: 



/tot — /tGB "h /disl + /ch 

Dd 2 



2lbld 
E 



n=l \ / a n _ 1 



+ 



hd 



hid hid ' 



(5.2) 



(5.3) 



where E is the energy cost per unit line of an i ndiv idual 
screw dislocation. In the free energy density ([5.2|), the 
total energy cost of dislocations given by the first two 
terms competes with the gain in the chiral energy. The 
twist penetratio n de pth A sets the length scale for Id and 
lb- Inspecting (p 7 ^), we see that the optimal values of 
ld/\ h/ A arc controlled by a single combination of the 
material parameters 



We minimized of the energy density (5.2) with respect 
to Zd/A and Z^/A numerically. The results are presented 
m Fig. | We can make several observations regarding 
our results. As is directly evident from Fig. |^, there is 
a range of a where the preferred values of Id and lb are 
of several A and, moreover, the ratio ld/h is close to 1. 
According to the remarks at the beginning of this sec- 
tion, this kind of geometry validates our computational 
techniques. To see how much faraway grain boundaries 
contribute to the establishment of the dislocation lattice 
structure, we computed the posi tion s of minima for free 
energy densities obtained from ( |5.2| ) by truncating the 
sums over n to one, two, and three terms. The results of 
this computation are presented in Fig. ^, which demon- 
strate that in the range < a < 0.2, which corresponds 
to lb, Id > 2.5A, the lattice structure is determined almost 







entirely by the interactions of nearest and next-to-nearest 
twist grain boundaries. 



a 



2(hd - E) 



(5.4) 



3.5 



2.5 




0.1 0.2 0.3 

a = 2(hd-E)/Dd 2 



0.4 



FIG. 3. Dependence of the minima positions on the num- 
ber of terms kept in the sum over grain boundaries in the 
free energy density (|3.2|). The style of each curve reflects the 
number of terms used in its computation. The upper set of 
curves correspond to ?<j/A, while the lower curves correspond 
to lb/\. The solid curves are the same curves as in Fig. ^. 



constructed by analogy with the theory of vortex inter- 
action in the Abrikosov phase in the high-K limit. The 
resulting theory was applied to the calculation of the lat- 
tice parameters of the screw dislocation arrangement. We 
found that in this limit the ratio of the two lattice param- 
eters lb /Id remains pratically constant over a very wide 
range of the control parameter. The value of the ratio is 
0.95. It is interesting to note this value is within several 
thousands of the value obtained by Renn and Lubensky 
jjj in the opposite limit h — > h C 2 for a specific value 
K/K2 = 0. Whether this is accidental or not remains to 
be discovered. In contrast to the Renn and Lubensky re- 
sult, the value oild/h obtained here is independent of Ki 
and K3. This value is consistent with the experimental 
data of Navaillcs and coworkers [fL5| . 
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0.15 0.2 0.25 0.3 



1.1 




0.9 



0.8 



0.96 , i 

0.958 
0.956 
~ 0.954 
0.952 

0.95 ' 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 

a = 2(hd-E)/Dd 2 

FIG. 4. Dependence of h/ld on the control parameter a. 
Note that it is nearly independent of a. 

A very interesting result emerges regarding the ratio of 
the lattice parameters lb /id- It turns out that the ratio 
is nearly constant for a very wide range of values of the 
control parameter a. This prediction should be compared 
to the experimental values of 4 /Id measured by Navailles 
and coworkers [Jl5||. 



VI. CONCLUSIONS 

The theory of the screw dislocation interaction in the 
TGB^ phase in the low-angle grain boundary limit was 



0.7 



100 100.5 101 101.5 

Temperature in °C 



102 



FIG. 5. Experimental dependence of h/ld on temperature. 
The data set marked with triangles is taken from |L5| , the data 
set marked with boxes was provided courtesy of L. Navailles. 
The first set of data were taken while increasing temperature, 
while the latter was taken in a run with decreasing tempera- 
ture. Note that the grain rotation angle increases from 6° for 
the lowest temperature to 9° for the highest temperature. 



Recent work |T^] which studied the nonlinear elastic- 
ity of smectic liquid crystals showed that defects in the 
same grain boundary had power-law interactions. One 
might expect, in general, that the interaction between 
grain boundaries would remain exponential. It would 
seem then that the grain boundaries would move closer 
together as the defects in the boundaries would move 
further apart, casting some doubt on our result (and on 
experiment). In the absence of director modes, the non- 
linear elasticity is [B: 



d 3 x 



B 



d z u - ^ (Vw) 2 



-i 2 



(6.1) 



In this nonlinear theory the director and the layer normal 
are locked together so that n = — V</>/|V0|. While a 
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full analysis of the energetics of two interacting grain 
boundaries will be the focus of further work, we can, in 



the spirit of the analysis at the end of section IV, consider 
a single, distorted screw dislocation: 



u(x, y) = tan 



Try 



I t&n(nx/l) 



(6.2) 



The layer normal of this dislocation relaxes in the usual 
way along the y-direction but relaxes to its asymptotic 
value in the confined region between x — ±1/2. It is 
straightforward to calculate the nonlinear energy of this 
defect and to find the "confinement" energy by subtract- 
ing the / = oo value. Expanding in powers of l^ 1 we 
have 



V 2 u 



tt 2 2xy(x 2 — y 2 ) 



V2 ( x 2 + y 2)2 



(Vti) 



■ir 



9 2 2 2 4 

tt x y — |cc 
P (x 2 + y 2 ) 2 ' 



(6.3) 
(6.4) 



It is clear that the ^-dependence of the confinement en- 
ergy scales as l~ 2 which can balance the l^ 2 interaction 
found between defects in the same grain boundary [ jT9"| . 
We shall explore this further in future work. 
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